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This  paper  presents  a  five-year  global  simulation  of  HYCOM.  the  HYbrid  Coordinate  Ocean  Model,  that 
simultaneously  resolves  the  eddying  general  circulation,  barotropic  tides,  and  baroclinic  tides  with  32 
layers  in  the  vertical  direction  and  1/12.5®  (equatorial)  horizontal  grid  spacing.  A  parameterized  topo¬ 
graphic  wave  drag  is  inserted  into  the  model  and  tuned  so  that  the  surface  tidal  elevations  are  of  com¬ 
parable  accuracy  to  those  in  optimally  tuned  forward  tide  models  used  in  previous  studies.  The  model 
captures  93%  of  the  open-ocean  sea-surface  height  variance  of  the  eight  largest  tidal  constituents,  as 
recorded  by  a  standard  set  of  102  pelagic  tide  gauges  spread  around  the  World  Ocean.  In  order  to  min¬ 
imize  the  impact  of  the  wave  drag  on  non-tidal  motions,  the  model  utilizes  a  running  25-h  average  to 
approximately  separate  tidal  and  non-tidal  components  of  the  near-bottom  flow.  In  contrast  to  earlier 
high-resolution  global  baroclinic  tide  simulations,  which  utilized  tidal  forcing  only,  the  simulation  pre¬ 
sented  here  has  a  horizontally  non-uniform  stratification,  supported  by  the  wind-  and  buoyancy  forcing. 
The  horizontally  varying  stratification  affects  the  baroclinic  tides  in  high  latitudes  to  first  order.  The  mag¬ 
nitude  of  the  internal  tide  perturbations  to  sea  surface  elevation  amplitude  and  phase  in  a  large  box  sur¬ 
rounding  Hawai'i  is  quite  similar  to  that  observed  in  satellite  altimeter  data,  although  the  exact  locations 
of  peaks  and  troughs  in  the  modeled  perturbations  differ  from  those  in  the  observed  perturbations. 

©  2010  Elsevier  Ltd.  All  rights  reserved. 


1.  Introduction 

This  paper  presents  an  early  attempt  to  simultaneously  resolve 
the  oceanic  general  circulation,  its  associated  mesoscale  eddy  field, 
and  the  barotropic  and  baroclinic  tides,  at  high  horizontal  and  ver¬ 
tical  resolution,  in  a  global  model.  A  number  of  scientific  and  oper¬ 
ational  topics  can  be  addressed  in  such  a  model-for  instance,  the 
interactions  between  tides  and  low-frequency  motions,  and  the  ef¬ 
fect  of  eddies  on  the  predictability  of  internal  tides.  We  forgo  some 
of  these  interesting  scientific  questions  here,  until  more  detailed 
and  time-consuming  analyses  of  our  results  are  undertaken.  The 
focus  here  is  on  the  numerical  techniques  we  have  used  to  ensure 
an  accurate  barotropic  tide  without  severely  disrupting  the  eddy¬ 
ing  general  circulation.  As  will  be  described  below,  it  is  far  from 
trivial  to  ensure  an  accurate  barotropic  tide  in  forward  global  mod¬ 
els.  and  the  presence  of  non-tidal  motions  only  increases  the  chal¬ 
lenge.  We  also  present  a  few  preliminary  results  from  our 
simulations.  We  present  comparisons  with  observations  from  both 
satellite  altimeters  and  tide  gauges.  These  comparisons  indicate 
that  we  have  achieved  some  success  in  reproducing  the  observed 
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barotropic  and  baroclinic  tide  fields.  We  demonstrate  the  first-or¬ 
der  effect  of  horizontally  varying  stratification  on  the  global 
internal  tides.  Finally,  we  provide  a  few  visual  demonstrations  of 
the  co-existence  of  barotropic  tides,  baroclinic  tides,  the  general 
circulation,  and  mesoscale  eddies  in  the  model. 

In  recent  years,  several  groups  have  simulated  the  global  oce¬ 
anic  general  circulation  in  numerical  models  with  horizontal  grids 
that  are  fine  enough  to  resolve  (or  at  least  permit)  mesoscale  ed¬ 
dies,  the  transient  turbulent  features  which  contain  a  substantial 
fraction  of  the  oceanic  kinetic  energy.  For  instance,  the  Parallel 
Ocean  Program  (POP)  model  has  been  run  globally  at  1/10°  resolu¬ 
tion  (Maltrud  and  McClean,  2005),  the  Naval  Research  Laboratory 
Layered  Ocean  Model  (NLOM)  is  run  in  ocean  forecast  mode  with 
1/32°  horizontal  resolution  (Shriver  et  al..  2007),  the  HYbrid  Coor¬ 
dinate  Ocean  Model  (HYCOM)  is  being  developed  as  a  1/12.5°  res¬ 
olution  forecast  model  (Chassignet  et  al.,  2007).  the  Ocean  General 
Circulation  Model  for  the  Earth  Simulator  (OfES)  has  been  run  at  1  / 
10°  resolution  (Masumoto  et  al..  2004).  and  the  Ocean  Circulation 
and  Climate  Advanced  Model  (OCCAM)  has  achieved  1/12°  hori¬ 
zontal  resolution  (Lee  et  al.,  2007).  At  the  same  time,  in  recent 
years,  high-resolution  global  models  of  the  baroclinic  tides  have 
begun  to  be  run  (Arbic  et  al..  2004  -  hereafter,  ACHS;  Simmons 
et  al..  2004  -  hereafter,  SHA:  Hibiya  et  al.,  2006;  Simmons.  2008). 
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In  coastal  models,  it  is  common  to  model  tides  and  non-tidal  mo¬ 
tions  simultaneously.  However,  tides  and  non-tidal  motions  have 
almost  always  been  simulated  separately  in  global  models.  A  few 
recent  global  simulations  have  included  tides  and  non-tidal  mo¬ 
tions  simultaneously  (Schiller  and  Fiedler,  2007;  T.  Dobslaw,  M. 
Muller,  and  M.  Thomas,  personal  communication.  2008),  but  these 
studies  are  done  with  model  horizontal  grid  spacings  of  order  one 
degree,  at  which  neither  mesoscale  eddies  nor  baroclinic  tides  are 
resolved.*  Here  we  merge  two  previously  separate  recent  threads  in 
the  literature  -  high-resolution  modelling  of  the  global  eddying  gen¬ 
eral  circulation,  and  high-resolution  modelling  of  the  global 
barotropic  and  baroclinic  tides. 

By  combining  these  two  threads  we  potentially  improve  the 
modelling  of  both  types  of  motions,  which  affect  each  other  in  var¬ 
ious  ways.  Interactions  between  mesoscale  eddies  and  internal 
tides  have  the  potential  to  transfer  part  of  the  coherent  internal 
tide  energy  into  incoherent  signals,  and  to  affect  tidal  energy  bud¬ 
gets  (Park  and  Watts.  2006;  Rainville  and  Pinkel,  2006:  Zaron  et  al., 
2009;  Chavanne  et  al.,  in  press-a,  in  press-b).  Park  and  Watts 
(2006)  and  Chavanne  et  al.  (in  press-b)  show  that  the  variations 
in  stratification  induced  by  mesoscale  eddies,  in  addition  to  the 
scattering  arising  from  eddy  velocities,  have  important  effects  on 
internal  tide  propagation.  A  mixed  tidal/non-tidal  model  is  also 
more  likely  to  properly  account  for  the  effects  of  the  quadratic  bot¬ 
tom  boundary  layer  drag  term.  Currently,  many  ocean  general  cir¬ 
culation  models  insert  an  assumed  tidal  background  flow,  typically 
taken  to  be  about  5cms"\  into  the  quadratic  drag  formulation 
(e.g.  Willebrand  et  al.,  2001).  However,  in  the  actual  ocean  tidal 
velocities  vary  from  -^1-2  cm  s"*  in  the  abyss,  to  ~50-100cm 
s  *  in  areas  of  large  coastal  tides.  Thus  an  assumed  tidal  back¬ 
ground  flow  of  5  cm  S'*  is  too  strong  in  the  abyss,  and  too  weak 
in  coastal  areas.  By  actually  resolving  the  (spatially  inhomoge¬ 
neous)  tidal  flows  in  a  general  circulation  model,  we  take  a  step  to¬ 
wards  correcting  this  problem.  The  explicit  resolution  of  tides  may 
represent  an  important  step  towards  more  realistic  representation 
of  mixing  in  high-resolution  models,  and  we  are  currently  pursuing 
this  avenue  as  well.  Finally,  the  stratification  in  a  mixed  tidal/non- 
tidal  model  can  vary  horizontally  as  in  the  actual  ocean,  since  the 
wind-  and  buoyancy-forcing  which  supports  this  varying  stratifica¬ 
tion  is  present.  In  contrast,  the  stratification  in  the  earlier  high-res- 
olution  global  baroclinic  tide  simulations  of  AGHS  and  SHA  was 
chosen  to  be  horizontally  uniform  since  these  simulations  did 
not  include  wind-  and  buoyancy-forcing. 

The  results  presented  here  represent  an  important  first  step  to¬ 
wards  one  of  our  long-term  goals,  to  simultaneously  resolve  tides 
and  non-tidal  motions  in  global  data-assimilative  models  with  1/ 
25°  horizontal  resolution.  Because  the  goal  is  an  operational  model, 
accuracy  of  all  the  resolved  motions  is  paramount.  We  therefore,  de¬ 
sire  to  begin  with  forward  tide  models  that  are  as  accurate  as  possi¬ 
ble.  In  recent  years,  it  has  been  shown  that  achieving  accurate 
surface  elevations  in  forward  global  barotropic  tide  models  requires 
the  insertion  of  a  parameterization  of  drag  (and  energy  loss)  due  to 
the  breaking  of  internal  waves  generated  by  tidal  flow  over  rough 
topography  (Jayne  and  St.  Laurent,  2001;  Carrere  and  Lyard,  2003; 
Egbert  etal., 2004;  AGHS;  Lyard  etal.,  2006;  Uehara  etal..2006;  Grif¬ 
fiths  and  Peltier,  2008,  2009).  These  parameterizations  are  moti¬ 
vated  by  inferences  from  tide  models  constrained  by  satellite 
altimetry  of  the  dissipation  of  tidal  energy  in  mid-ocean  areas  of 
rough  topography  (Egbert  and  Ray.  2000),  as  well  as  in-situ  evidence 
of  elevated  dissipation  levels  in  such  areas  (e.g.  Polzin  et  al.,  1997). 
The  subtleties  of  applying  a  parameterized  topographic  wave  drag 
in  models  which  resolve  the  generation  of  low-mode  baroclinic 


'  To  be  more  precise  about  the  Schiller  and  Fiedler  (2007)  simulation,  their 
resolution  was  high  in  an  area  around  Australia,  but  the  telescoping  grid  they  used  led 
to  low  resolutions  over  most  of  the  global  ocean. 


tides,  and  in  models  which  resolve  non-tidal  as  well  as  tidal  motions, 
will  be  discussed  in  the  next  section.  A  comparison  and  discussion  of 
the  accuracies  of  the  barotropic  tides  in  the  baroclinic  simulations  of 
AGHS  and  SHA  will  prove  to  be  instructive  with  regard  to  handling 
topographic  wave  drag  in  the  new  HYCOM  simulations. 


2.  Inclusion  of  parameterized  topographic  wave  drag 

2.1.  Need /or  poTameterized  wave  drag  in  barac/inic  ride  models 

In  barotropic  tide  models,  none  of  the  internal  waves  generated 
by  flow  over  rough  topography  are  resolved,  and  all  of  this  wave 
activity  must  be  parameterized.  In  baroclinic  tide  models,  the  situ¬ 
ation  is  more  complicated  and  interesting.  The  resolved  generation 
of  low-mode  internal  tides  means  that  the  barotropic  tide  will  be 
losing  energy  to  the  baroclinic  tide  in  baroclinic  models.  Indeed, 
the  computation  of  this  energy  conversion  was  a  central  goal  of 
SHA,  which  built  upon  the  baroclinic  tide  simulations  performed 
for  AGHS.  Both  studies  were  done  with  HIM.  the  Hallberg  Isopycnal 
Model  (Hallberg  and  Rhines,  1996).  Since  in  baroclinic  tide  models 
energy  is  lost  from  the  barotropic  mode,  it  is  tempting  to  view 
parameterized  topographic  wave  drag  as  redundant.  Indeed,  the 
main  baroclinic  simulation  of  SHA  did  not  retain  the  parameterized 
topographic  wave  drag  used  in  the  main  AGHS  baroclinic  simula¬ 
tions.  We  now  examine  the  consequences  of  these  different  choices 
made  in  AGHS  and  SHA. 

Table  1  shows  the  globally  integrated  available  potential  en¬ 
ergy  (APE)  at  the  sea  surface,  and  the  globally  integrated  baro¬ 
tropic  kinetic  energy  (KE).  both  computed  via  standard  formulae 
which  can  be  looked  up  in  for  instance  AGHS.  in  (1)  the  satellite 
altimeter-constrained  barotropic  solutions  of  Egbert  and  Ray 
(2003),  (2)  the  main  baroclinic  simulation  of  AGHS  (see  their 
Fig.  11),  which  utilized  parameterized  topographic  wave  drag 
optimally  tuned  to  minimize  sea  surface  elevation  errors  with  re¬ 
spect  to  satellite  altimetry.  (3)  a  baroclinic  simulation  of  AGHS 
which  did  not  utilize  parameterized  topographic  wave  drag  and 
which  also  used  only  the  scalar  approximation  (e.g.  Ray.  1998) 
for  the  self-attraction  and  loading  term  (in  other  words,  run  un¬ 
der  conditions  similar  to  the  main  baroclinic  simulation  of 
SHA),  (4)  the  main  baroclinic  simulation  of  SHA,  and  (5)  a  baro¬ 
clinic  simulation  of  SHA  briefly  mentioned  in  their  appendix,  in 
which,  inspired  by  Fig.  2  of  AGHS.  an  unrealistically  large  value 
of  Cd  (100  times  the  normal  value)  was  utilized  as  a  proxy  for 
topographic  wave  drag.  The  globally-  and  temporally-averaged 
root-mean-square  (RMS)  elevation  errors  of  the  forward  models 
with  respect  to  GOT99  (Ray.  1999),  a  highly  accurate  altimetry- 
constrained  tide  model,  are  also  shown.  The  errors  are  computed 
over  waters  deeper  than  1000  m  and  over  latitudes  covered  by 
the  TOPEX/POSEIDON  altimeter  (equatorward  of  66°).^  Finally, 
the  percentage  of  the  GOT99  open-ocean  sea  surface  height  (SSH) 
variance  captured  by  the  models  is  shown.  AGHS  may  be  consulted 
for  details  of  how  the  errors  and  percent  variance  captured  are  cal¬ 
culated.  In  the  main  baroclinic  simulation  of  AGHS,  the  surface  APE 
and  barotropic  KE  are  both  quite  close  to  the  Egbert  and  Ray  val¬ 
ues,  and  as  a  result  the  surface  elevation  error  is  reasonably  small. 
Note  also  that  the  barotropic  energies  and  elevation  errors  in  this 
optimally  tuned  AGHS  baroclinic  simulation  are  barely  different 
from  those  in  the  optimally  tuned  one-layer  simulation  of  AGHS 
(not  shown).  On  the  other  hand,  in  the  AGHS  baroclinic  simulation 
run  without  any  parameterized  topographic  wave  drag,  the  surface 
APE  and  barotropic  KE  are  both  about  twice  as  large  as  the  ob- 


*  To  be  precise,  the  errors  in  the  SHA  results  were  computed  over  the  btitude  range 
66*^  to  64'’N.  in  order  to  avoid  the  complex  tripolar  grid  utilized  in  the  high  latitudes 
of  that  study  in  the  error  computations. 
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Energies  and  elevation  errors  of  earlier  forward  global  baroclinic  tide  models, 
compared  to  results  from  an  altimetry-constrained  tide  model.  Globally  integrated 
surface  available  potential  energy  (APE)  and  barotropic  kinetic  energy  (KE)  of  Mj  are 
computed  from  ( I )  Table  I  of  Egbert  and  Ray  (2003  -  ER2003  below-based  on  an 
altimetry-constrained  barotropic  model),  (2)  main  baroclinic  simulation  of  AGHS  (see 
their  Fig.  11).  (3)  AGHS  baroclinic  simulation  without  parameterized  topographic 
wave  drag  and  with  only  a  scalar  approximation  for  self-attraction  and  loading  (SAL); 
i.e.  conditions  like  those  in  the  main  baroclinic  simulation  of  SHA,  (4)  main  SHA 
baroclinic  simulation  (see  their  Fig.  8),  (5)  a  baroclinic  simulation  in  the  appendix  of 
SHA.  with  an  artificially  large  c^  value  of  0.3  standing  in  as  a  proxy  for  parameterized 
topographic  wave  drag.  Units  of  energies  are  lO”  J.  Globally  averaged  sea-surface 
elevation  discrepancies  (cm)  with  respect  to  GOT99,  computed  over  latitudes 
equatorward  of  66"  and  waters  deeper  than  1000  m,  are  also  given.  Numbers  in 
parentheses  indicate  percent  of  altimeter-measured  open-ocean  sea  surface  elevation 
variance  captured  by  the  models. 


Model 

Surface 

APE 

Barotropic 

KE 

RMS  elevation 

error 

ER2003 

134 

1.78 

- 

AGHS  main  baroclinic 
simulation 

1.48 

1.73 

7.37  (92.4) 

AGHS.  no  parameterized 
drag,  scalar  SAL 

3.18 

3.46 

17.14(58.8) 

SHA  main  baroclinic 
simulation 

4.37 

5.09 

23.35  (23.5) 

SHA.  large  c^  as  proxy 
for  wave  drag 

1.66 

2.03 

9.88  (86.3) 

served  values.  As  a  consequence  the  elevation  discrepancy  with  re¬ 
spect  to  GOT99  is  much  larger,  and  the  percentage  of  SSH  variance 
captured  is  much  lower.  Consistent  with  this  result,  both  the  po¬ 
tential  and  kinetic  energies  of  the  main  SHA  baroclinic  simulation 
are  also  larger,  by  factors  of  about  3,  than  those  in  the  accurate  sa¬ 
tellite-constrained  models,  and  the  high  elevation  error  and  low 
percent  variance  captured  reflect  this  mismatch.^  The  large  c^  sim¬ 
ulation  mentioned  in  the  appendix  of  SHA  performs  much  better 
with  respect  to  the  observations  (and  sees  a  factor  of  2.4  drop  in 
the  conversion  of  barotropic  to  low-mode  baroclinic  energy),  dem¬ 
onstrating  that  even  artificial  frictions  can  lead  to  accurate  modeled 
barotropic  tides  as  long  as  they  remove  energy  at  approximately 
the  correct  rate. 

In  the  current  study  we  have  also  found  that  the  barotropic 
tide  is  extremely  inaccurate  if  parameterized  topographic  wave 
drag  is  not  included.  As  in  AGHS,  we  choose  to  have  the  topo¬ 
graphic  wave  drag  in  HYCOM  acting  on  the  near-bottom  flow 
(in  the  HYCOM  multi-layer  simulations  we  define  ’’near-bottom" 
flow  as  flow  averaged  over  the  bottom  500  m).  AGHS  argued  that 
a  conversion  of  energy  from  barotropic  to  baroclinic  tides  in 
baroclinic  models  does  not  represent  a  loss  of  energy  in  the  total 
(barotropic  plus  baroclinic)  system,  since  the  models  do  not  re¬ 
solve  the  breaking  of  baroclinic  tides  occuring  in  the  actual 
ocean.  Parameterized  topographic  wave  drag  acting  on  the 
near-bottom  flow,  unlike  modal  conversion,  drains  total  energy 
from  the  model.  The  near-bottom  flow  is  a  function  of  both 
barotropic  and  baroclinic  tides,  the  latter  contributing  less  when 
the  stratification  is  surface-intensified,  as  is  typical  in  the  ocean. 
If  we  assume  that  wave  breaking  takes  place  mostly  in  the  deep 
ocean,  just  above  rough  topography,  and  involves  mostly  high 
vertical  modes,  then  the  parameterization  represents  the  break¬ 
ing  of  high  modes  near  the  bottom,  which  is  not  resolved  in 
present-day  global  baroclinic  tide  models.  It  remains  to  be  seen 
whether  this  is  the  best  representation  of  what  actually  happens 
in  nature  (some  of  the  tidal  energy  loss  in  the  actual  ocean  may 


’  Consistent  with  these  results  of  AGHS  and  SHA,  a  new  simulation  of  baroclinic 
tides  performed  with  ^50  z-levels  in  the  vertical  direction  and  1/4"  horizontal 
resolution  shows  excessively  large  barotropic  tides  in  the  absence  of  parameterized 
topographic  wave  drag  (Andrew  Coward.  Ariane  Koch-Larrouy.  Gurvan  Medec,  Adrian 
New.  George  Nurser,  and  David  Smeed.  personal  communication.  2009). 


be  due  to  high  vertical  mode  breaking  in  the  upper  ocean),  but  it 
is  clearly  true  empirically  that  it  results  in  far  superior  barotrop¬ 
ic  tides  than  those  in  simulations  which  do  not  utilize  parame¬ 
terized  topographic  wave  drag.  It  is  clear  that  any  forward 
global  tide  model  that  aspires  to  be  the  starting  point  for  an 
operational  model  must  include  a  parameterized  topographic 
wave  drag,  or  some  other  way  of  removing  total  energy  from 
the  model. 

The  need  for  a  parameterized  topographic  wave  drag  to  achieve 
accurate  barotropic  tides  is  linked  with  the  horizontal  and  vertical 
resolutions  available  in  the  model;  if  the  breaking  of  high  vertical 
mode  waves  and  consequent  energy  loss  could  be  properly  re¬ 
solved,  then  parameterized  topographic  wave  drag  would  not  be 
needed.  The  HYCOM  simulations  presented  here  utilize  horizontal 
and  vertical  resolutions  that  are  state-of-the-art,  in  accordance 
with  the  largest  computers  presently  available.  Prompted  by  the 
suggestion  of  one  of  the  reviewers  (Sonya  Legg),  we  briefly  discuss 
here  the  limitations  of  these  resolutions.  From  Fig.  11  of  SHA.  we 
take  the  wavelength  of  the  first  baroclinic  mode  M2  internal  tide 
to  be  about  100  km  in  mid-latitudes.  A  widely  used  rule  of  thumb 
holds  that  eight  gridpoints  per  wavelength  are  required  to  properly 
resolve  waves  in  a  model.  The  1/12.5°  horizontal  resolution  in  our 
simulations  translates  to  about  9  km  grid  spacing.  Therefore,  the 
mode  one  internal  tides  are  well-resolved  in  our  simulations.  How¬ 
ever,  since  the  wavelengths  drop  off  approximately  as  1  /n.  where  n 
is  the  vertical  mode  number,  then  the  mode  two  internal  wave  has  a 
wavelength  of  about  50  km.  meaning  there  are  less  than  8  grid- 
points  per  wavelength.  By  the  time  we  reach  vertical  mode  num¬ 
bers  of  about  ten  or  so,  there  will  be  only  one  gridpK)int  per 
wavelength.  All  of  this  suggests  that  vertical  mode  numbers  beyond 
about  10  are  probably  not  resolved  at  all  in  the  simulations  pre¬ 
sented  here,  and  vertical  mode  numbers  beyond  one  or  two  are 
probably  not  well-resolved.  Thus  horizontal  resolution  limitations 
are  in  part  responsible  for  the  fact  that  parameterized  topographic 
wave  drag  is  still  required  to  achieve  accurate  barotropic  tides  in 
baroclinic  tide  models. 

To  end  this  subsection  we  note  that  despite  the  fact  that  the 
low-mode  baroclinic  tides  have  a  weak  signature  at  the  bottom, 
it  is  evident  that  insertion  of  parameterized  topographic  wave  drag 
into  a  baroclinic  tide  model  affects  the  propagation  distances  of  the 
low-mode  internal  tides.  Contrast,  for  instance,  the  shorter  propa¬ 
gation  distances  of  the  low-mode  internal  tide  beams  from  their 
source  regions  shown  in  Fig.  1 1  of  AGHS  with  the  longer  distances 
seen  in  Fig.  8  of  SHA. 


2.2.  Adapfafion  0/  paramererized  wave  drag  used  in  previous  sfudies 

We  utilize  an  adaptation  of  the  topographic  wave  drag  scheme 
described  in  the  appendix  of  AGHS.  which  is  based  on  the  scheme 
outlined  in  Garner  (2005).  A  multiplicative  factor  was  included  in 
the  scheme  and  tuned  to  minimize  the  globally  averaged  deep- 
ocean  RMS  elevation  discrepancy  between  the  forward  model 
and  GOT99.  AGHS  found  that  the  optimal  multiplicative  factor, 
for  1/2°  simulations  with  wave  drag  acting  only  in  waters  deeper 
than  1000  m,  was  about  7.  AGHS  suggested  that  the  multiplicative 
factor  may  compensate  for  the  small  scales  that  are  absent  in  the 
roughness  of  present-day  topographic  datasets  (e.g.  Smith  and 
Sandwell,  1997).  This  suggestion  will  shortly  be  tested,  in  tide 
model  runs  which  utilize  global  grids  of  statistical  roughness  pro¬ 
duced  by  Goff  and  Arbic  (2010).  These  grids  utilize  relationships 
between  geophysical  properties  such  as  seafloor  spreading  rate 
and  abyssal  hill  statistics  to  produce  synthetic  topographies.  The 


*  Along  with  several  collaborators,  we  are  exploring  the  possibility  of  removing 
energy  directly  from  the  resolved  vertical  shear. 
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Fig.  1.  Enfolding  time  (days)  for  parameterized  topographic  wave  drag  with  drgscl  -  1.  The  majority  of  the  white  areas  have  no  wave  drag  at  all. 


statistical  roughness  grids  can  be  overlaid  on  the  Smith  and  San- 
dwell  (1997)  dataset  to  create  global  bathymetries  with  the  right 
statistical  texture  of  the  seafloor,  though  they  are  not  determinis¬ 
tically  correct.  We  can  use  the  new  rougher  topography  as  an  input 
to  our  drag  scheme  to  determine  if  smaller  multiplicative  factors 
can  then  be  used.  It  is  important  to  note  that  while  these  exercises 
will  test  the  robustness  of  the  underlying  wave  drag  scheme  used 
in  AGHS  (since  a  multiplicative  factor  of  7  should  not  be  necessary 
for  a  correct  scheme  that  utilizes  a  topography  of  adequate  resolu¬ 
tion).  they  will  not  greatly  affect  the  total  drag  required  by  the 
model  to  bring  about  accurate  tides.  As  noted  by  AGHS,  the  glob¬ 
ally  averaged  strength  of  wave  drag  used  in  AGHS  and  in  Jayne 
and  St.  Laurent  (2001),  which  utilized  completely  different 
schemes  and  tunable  parameters,  was  virtually  identical.  As  a  prac¬ 
tical  matter,  a  certain  amount  of  drag  is  needed  to  obtain  accurate 
barotropic  tides.  Whether  the  underlying  theory  can  come  to  this 
correct  amount  of  drag  without  a  tunable  parameter  such  as  our 
multiplicative  factor  is  an  important  question,  but  one  of  limited 
practical  consequence. 

For  the  sake  of  simplicity,  here  we  reduce  the  AGHS  tensor 
scheme  to  a  scalar  scheme,  utilizing  energy  considerations.  We 
compute  from  1/8®  runs  of  the  AGHS  model  the  quantity 

^  Iropodrqg  ' 

”  (U‘U)  '  ^ 

where  angle  brackets  denote  time-averaging,  u  is  the  velocity  vec¬ 
tor,  and  the  term  in  the  momentum  equation  of  the 

(barotropic)  AGHS  model  arising  from  the  full  tensor  form  of  the 
topographic  wave  drag.  We  use  the  resulting  map  of  r  in  the  HYCOM 
simulations.  Note  that  r  is  a  linear  drag  coefficient,  with  1  /r  as  its  e- 
folding  time  scale.  Here  we  set  the  values  of  r  to  zero  in  regions 
shallower  than  500  m  and  where  it  is  small,  i.e.  where  the  e-folding 
time  is  greater  than  30  days,  which  together  account  for  73%  of  the 
area  of  the  World  Ocean.  We  also  clip  the  value  of  r  so  that  its  min¬ 
imum  e-folding  time  is  9  h.  Both  of  these  actions  limit  the  impact  of 
the  wave  drag  on  non-tidal  motions  (see  next  section).  Fig.  1  shows 
maps  of  1  /r  values  obtained  after  all  of  these  changes  have  been 
implemented.  The  drag  is  concentrated  over  well-known  areas  of 


rough  topography  such  as  the  Mid-Atlantic  Ridge,  Southwest  Indian 
Ridge,  etc.  In  the  1/12.5®  simulations  we  will  be  presenting  shortly, 
we  found  that  a  multiplicative  factor  of  6  yielded  accurate  tides, 
similar  to  the  optimal  multiplicative  factor  of  about  7  found  by 
AGHS.^  Hence,  the  e-folding  time  of  the  applied  drag  is  1.5  h  to  six 
days  with  no  drag  over  73%  of  the  worlds  oceans. 

2.3.  Utilizing  fopograp/iic  wave  drag  in  the  presence  of  non-tidal 
motions 

On  the  relatively  fast  timescales  of  internal  gravity  waves,  low- 
frequency  motions  such  as  mesoscale  eddies  and  strong  currents 
(i.e.  the  Gulf  Stream,  Antarctic  Circumpolar  Current,  etc.)  can  be  re¬ 
garded  as  steady.  The  generation  of  internal  gravity  waves  (lee 
waves)  by  steady  flows  over  rough  topography  is  a  classic  problem 
in  geophysical  fluid  dynamics  (e.g.  Gill,  1982).  Tidal  motions  are 
oscillatory,  not  steady,  and  the  work  of  Bell  (1975)  shows  that 
the  wave  drag  resulting  from  oscillatory  flow  over  rough  topogra¬ 
phy  differs  from  the  wave  drag  resulting  from  steady  flow.  In  the 
future  we  may  wish  to  include  a  parameterized  wave  drag  for 
the  non-tidal  (steady,  in  this  context)  flow  over  rough  topography 
in  HYCOM.  Indeed,  some  recent  papers  have  argued  that  this 
mechanism  represents  a  substantial  energy  loss  for  low-frequency 
motions  (Naveira-Garabato  et  al.,  200A\  Marshall  and  Naveira- 
Garabato,  2008;  Nikurashin.  2008).  For  now,  however,  we  wish  to 
have  the  wave  drag  acting  only  on  the  tidal  part  of  the  flow.  This 
presents  a  challenge:  how  is  the  model  to  know  the  partition  of  ti¬ 
dal  versus  non-tidal  near-bottom  flows?  In  order  to  accomplish 
this  separation,  at  least  roughly,  we  utilize  running  25-h  averages 
before  applying  the  topographic  wave  drag  to  the  near-bottom 
flow.  The  details  of  this  scheme  are  discussed  next. 


*  Arbic  et  ^1.  (2008)  found  that  when  the  HIM  tide  model  was  run  at  higher 
resolution  (1/8°)  and  the  topographic  wave  drag  acted  in  shallow  as  well  as  deep 
waters,  accurate  tides  could  be  obtained  with  a  multiplicative  factor  as  small  as  3.  The 
higher  multiplicative  factor  needed  in  high-resolution  HYCOM  may  be  due  In  part  to 
the  clipping  described  earlier,  and  to  the  fact  that  the  wave  drag  does  not  operate  in 
shallow  waters.  The  optimal  multiplicative  factor  apparently  varies  by  a  factor  of 
about  2,  depending  on  model  resolution,  cutoff  depth  for  utilizing  the  wave  drag,  and 
other  factors. 
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2.4.  Separation  of  tidal  from  non- tidal  near-bottom  flows  in  wave  drag 
scheme 

The  topographic  wave  drag  is  nominally  applied  to  the  tidal 
flow  only  and  acts  on  the  bottom  500  m  of  the  water  column  (recall 
that  the  wave  drag  is  zero  in  waters  shallower  than  500  m).  To  fil¬ 
ter  out  the  tides  we  first  form  the  average  of  25  hourly  samples  of 
the  velocity  over  the  bottom  500  m  lagged  in  time  (i.e.  from  the 
previous  25  h).  Since  the  25-h  period  is  very  close  to  twice  that 
of  the  dominant  tidal  constituent  Mj,  most  of  the  tidal  motions 
are  filtered  out  of  these  averages.  This  is  the  detided  bottom  flow 
U(,  which  is  used  as  a  correction  to  standard  implicit  linear  friction 
over  the  bottom  500  m  of  the  water  column: 

gf+At  At  _  2Af(T(u'^^  -  5^)  (2) 

Here,  <t  is  the  linear  drag  coefficient  which  in  this  case  is  6r,  The 
friction  is  implicit  for  stability,  given  the  large  drag  coefficient.  A 
25-h  average  is  not  an  exact  tidal  filter,  and  lagging  it  in  time 
may  introduce  aliasing,  but  there  are  limits  on  what  is  practical  in 
a  running  ocean  model.  In  an  8-constituent  tides  only  test  case,  add¬ 
ing  the  25-h  filter  had  minimal  effect  on  the  accuracy  of  the  tides. 
This  issue  is  further  explored  in  Section  5.3.  It  is  important  to  reit¬ 
erate  that,  in  the  governing  equations  of  the  model,  the  separation 
of  tidal  from  non-tidal  flows  with  a  25-h  filter  is  done  only  during 
the  application  of  the  topographic  wave  drag  to  the  bottom  flow.  Ti¬ 
dal  and  low-frequency  flows  can  still  interact  with  each  other 
through  various  terms  in  the  model  equations-advection,  nonlin¬ 
earities  in  the  mass  conservation  equation,  and  quadratic  bottom 
boundary  layer  drag,  for  example.  Tidal  and  non-tidal  flows  are 
not  separated  from  each  other  in  these  other  cases  in  our  governing 
model  equations. 

3.  Implementation  of  self-attraction  and  loading 

Hendershott  (1972)  showed  that  global  numerical  tide  models 
must  account  for  self-gravitation  of  the  ocean  tide,  solid  earth 
deformation  due  to  the  load  of  the  ocean  tide,  and  perturbations 
to  the  gravitational  potential  due  to  the  self-gravitation  of  the  solid 
earth  thus  deformed.  Collectively,  these  terms  are  known  as  the 
self-attraction  and  loading  (SAL)  term.  A  complete  treatment  of 
the  SAL  term  requires  computing  a  spherical  harmonic  decomposi¬ 
tion  of  the  ocean  tide.  This  is  not  computationally  feasible  to  do  in 
the  model  as  it  runs,  and  instead  is  often  done  offline.  An  iterative 
procedure  appears  to  be  necessary  to  achieve  numerical  conver¬ 
gence  (e.g.  Egbert  et  al.,  2004,  ACHS).  In  the  model  runs  presented 
here,  as  was  done  in  SHA,  we  use  the  scalar  approximation,  in 
which  the  SAL  term  is  approximated  as  a  constant  p  times  the 
sea  surface  elevation  field  rj.  The  scalar  approximation  is  less  accu¬ 
rate  than  the  proper  spherical  harmonic  treatment,  but  is  far  more 
computationally  expedient.  In  tests  with  HYCOM  forced  only  by 
tides,  we  found  that  the  optimal  value  of  this  constant  in  terms 
of  minimizing  the  globally  averaged  RMS  sea  surface  elevation  dis¬ 
crepancy  with  GOT99  is  0.06.  Since  the  load  numbers  used  in  prop¬ 
er  calculations  of  the  SAL  term  depend  on  the  degree  of  spherical 
harmonic,  this  value  of  p,  which  was  determined  for  the  barotropic 
tides,  is  not  valid  for  smaller  scale  features  such  as  mesoscale  ed¬ 
dies  and  internal  tides. 

As  pointed  out  by  Hendershott  (1972),  the  SAL  term  should  ap¬ 
ply  to  non-tidal  as  well  as  tidal  flows.  However,  as  noted  above,  the 
value  of  the  SAL  scalar  p  we  use  is  not  appropriate  for  mesoscale 
eddies.  We  choose  to  apply  the  scalar  SAL  approximation  to  the 
non-steric  SSH,  which  is  dominated  by  barotropic  tides  (see 
Fig.  8).  The  non-steric  height  is  the  difference  between  the  SSH 
and  the  steric  SSH.  The  steric  SSH  anomaly  is  taken  to  be  the 
change  in  SSH  due  to  the  difference  between  the  instantaneous 


vertical  potential  density  profile  and  a  multi-year  model  mean  pro¬ 
file,  and  the  total  steric  SSH  is  taken  to  be  this  anomaly  plus  the 
corresponding  multi-year  model  mean  SSH.  This  procedure  en¬ 
sures  that  SAL  is  applied  to  the  barotropic  tides  at  the  same  time 
that  it  is  not  improperly  applied  (i.e.  with  an  inappropriate  p  value) 
to  mesoscale  eddies.  As  with  topographic  wave  drag,  the  problem 
of  applying  SAL  to  mixed  tidal/general  circulation  runs  deserves 
careful  consideration,  and  we  will  continue  to  explore  different  op¬ 
tions  for  the  SAL  term  in  future  work. 


4.  Other  details  of  the  HYCOM  simulations 

HYCOM  is  a  community  ocean  model  (http://oceanmodel- 
ing.rsmas.miami.edu/hycom/)  and  uses  a  generalized  (hybrid  iso- 
pycnal/terrain  -  following  ((7)/z-level)  vertical  coordinate  (Bleck, 
2002).  Typically,  the  model  includes  isopycnal  coordinates  in  the 
stratified  ocean  but  uses  the  layered  continuity  equation  to  make 
a  dynamically  smooth  transition  to  z-levels  (fixed-depth  coordi¬ 
nates)  in  the  unstratified  surface  mixed  layer  or  to  (7-levels  (ter¬ 
rain-following  coordinates)  in  shallow  water.  The  optimal 
coordinate  is  chosen  every  time  step  using  a  hybrid  coordinate 
generator.  In  this  way,  the  model  automatically  generates  the 
lighter  isopycnal  layers  needed  for  the  pycnocline  during  summer, 
while  the  same  layers  may  define  z-levels  during  winter. 

In  simulations  done  for  this  paper  HYCOM  uses  the  K-Profile 
Parameterization  (KPP;  Large  et  al.,  1994)  for  full-column  mixing, 
so  interior  diapycnal  mixing  is  from  KPPs  background  mixing 
terms.  The  background/intemal  wave  viscosity  is  set  to  1  cm^  s“^ 
and  the  diffusivity  to  0.1  cm^  s~\ 

The  model  spans  the  entire  globe  north  of  78.6“S,  with  a  Merca¬ 
tor  grid  from  66" S  to  47°N,  at  a  resolution  of  0.08^  cos(lat) 
xO  08'’(latitude  x  longitude),  and  a  bipolar  Arctic  patch  north  of 
47'’N,  i.e.  a  tripole  grid  (Murray,  1996).  The  average  zonal  (longitu¬ 
dinal)  resolution  for  this  1/12.5*  global  grid  varies  from  jss  9  km  at 
the  equator  to  w  7  km  at  mid-latitudes  (e.g.  at  40’N)  and  «  3.5  km 
at  the  north  pole.  The  meridional  (latitudinal)  grid  resolution  is 
halved  in  the  Antarctic  for  computational  efficiency.  The  total 
number  of  gridpoints  for  each  vertical  level  in  the  simulation  is 
3297  (nominally,  the  north-south  direction)  by  4500  (nominally, 
the  east-west  direction).  The  model’s  land-sea  boundary  is  at  the 
10-m  isobath  and  it  potentially  uses  a  terrain-following  vertical 
coordinate  in  depths  shallower  than  140  m.  The  bottom  topogra¬ 
phy  was  constructed  from  the  NRL  Digital  Bathymetry  Data  Base 
(DBDB2),  which  has  a  resolution  of  2-min  and  is  available  online 
at  http.7/www7320.nrIssc.navy.mil/DBDB2_WWW/.  Numerous 
hand-edits  have  been  performed  to  improve  coastlines  and  sill 
depths  in  key  straits  and  passages. 

There  are  32  hybrid  layers  in  the  vertical  direction  in  the  model. 
The  target  density  values  for  the  isopycnals  and  the  decreasing 
change  in  density  with  depth  between  isopycnal  coordinate  sur¬ 
faces  are  based  on  the  1/4°  Generalized  Digital  Environmental 
Model  (GDEM)  climatology  (NAVOCEANO,  2003;  Carnes,  2009). 
We  take  the  GDEM  z-level  climatology  to  create  an  ‘'isopycnal"  cli¬ 
matology  to  initialize  HYCOM  based  on  the  target  densities.  The 
density  difference  values  were  chosen,  so  that  the  layers  tend  to 
become  thicker  with  increasing  depth,  with  the  lowest  abyssal 
layer  being  the  thickest.  The  layer  thicknesses  are  not  the  same 
everywhere,  but  vary  spatially  (see  Fig.  10).  Upon  initialization,  a 
given  gridpoint  may  be  in  any  of  the  three  coordinate  types  (iso¬ 
pycnal,  z-level,  or  terrain  following)  and  this  may  change  with 
depth  as  well.  The  hybrid  coordinate  generator  favors  isopycnal 
layers.  The  minimum  thickness  of  each  layer  is  not  zero,  as  it  would 
be  for  an  outcropped  layer  in  a  purely  isopycnal  model,  but  rather  a 
fixed  value  per  layer  that  implicitly  creates  the  z-level  coordinates. 
So,  isopycnal  layers  effectively  outcrop  into  the  z-levels  formed  by 
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the  layers  above  rather  than  at  the  surface.  The  minimum  thick¬ 
ness  of  the  top  layer  in  deep  water  is  3  m.  and  this  minimum  in¬ 
creases  1.18x  per  layer  up  to  a  maximum  of  450  m,  and  target 
densities  are  chosen  such  that  at  least  the  top  four  layers  are  al¬ 
ways  in  z-level  coordinates. 

The  initial  model  spin-up  run  was  initialized  from  the  January 
CDEM  climatology  and  forced  by  years  1979-2002  from  the  Euro¬ 
pean  Centre  for  Medium-Range  Weather  Forecasts  (ECMWF)  40- 
year  Re-Analyses  (ERA-40)  (K^llberg  et  al.,  2004)  averaged  to  form 
a  climatological  monthly  mean  atmospheric  forcing.  The  wind 
speeds  were  scaled  to  be  consistent  with  QuikSCAT  observations 
(Kara  et  al.,  2009).  The  6-hourly  sub-monthly  wind  anomalies  from 
the  0.5®  Fleet  Numerical  Meteorology  and  Oceanography  Center 
(FNMOC)  Navy  Operational  Global  Atmospheric  Prediction  System 
(NOCAPS;  Rosmond  et  al..  2002)  over  year  2003  were  added  to  the 
12  monthly  averages  to  obtain  realistic  mixed  layer  depths  and  to 
allow  continuation  with  3-hourly  or  6-hourly  interannual  winds 
data  sets. 

5.  Results 

5.1.  Description  of  simulations  and  sampling  issues 

The  results  here  are  taken  from  five  different  simulations  of  HY- 
COM,  designated  by  9.7,  13.1,  14.0,  14.1,  and  14.2  (Table  2).  HY- 
COM  13.1  (’Tideonly**)  is  a  simulation  with  a  horizontally 
uniform  two  layer  stratification  and  only  tidal  forcing,  performed 
so  that  comparison  with  wind-,  buoyancy-,  and  tidal  ly  forced  sim¬ 
ulations  will  reveal  effects  of  horizontally  non-uniform  stratifica¬ 
tion  on  internal  tides.  HYCOM  9.7  serves  as  our  ’’ControP 
experiment,  with  wind-  and  buoyancy-forcing  but  without  any  ti¬ 
dal  forcing.  It  started  from  the  end  of  the  spin-up  simulation  and 
was  run  from  2003  through  mid-2008  using  3-hourly  FNMOC  NO¬ 
CAPS  atmospheric  forcing  with  wind  speeds  scaled  to  be  consistent 
with  QuikSCAT  observations.  HYCOM  14.0  (’‘M2tide")  was  a  test 
experiment,  performed  for  just  two  months  starting  in  July  2003 
from  HYCOM  9.7  (Control).  HYCOM  14.0  (M2tide)  included  Mj  ti¬ 
dal  forcing  as  well  as  wind-  and  buoyancy-forcing.  Encouraged 
by  the  results  of  14.0  (M2tide),  we  then  proceeded  to  14.1  (“All- 
tides’*),  an  experiment  again  starting  from  9.7  (Control)  in  July 
2003  but  covering  5  calendar  years  after  allowance  for  spinup 
(2004-2008).  HYCOM  14.1  (Alltides)  included  tidal  forcing  for 
M2,  S2.  N2.  and  K2  (the  four  largest  semidiurnal  constituents),  and 

Table  2 


Simulation  number,  shorthand 
described  in  this  paper. 

name,  and  characteristics  of  HYCOM  simulations 

Simulation 

Shorthand 

Characteristics 

number 

name 

9.7 

Control 

32-layers 

Wind-  and  buoyancy  forcing 

No  tidal  forcing  or  topographic  wave  drag 

13.1 

Tideonly 

2-layers 

No  wind-  and  buoyancy  forcing 

M2  tidal  forcing  and  topographic  wave  drag 

14.0 

M2tide 

32-layers 

Wind-  and  buoyancy  forcing 

M2  tidal  forcing  and  topographic  wave  drag 
Short  duration  (two  months) 

14.1 

Alltides 

32-layers 

Wind-  and  buoyancy  forcing 

Eight-constituent  tidal  forcing  and 
topographic  wave  drag 

14.2 

3Doutput 

32 -layers 

Wind-  and  buoyancy  forcing 

Eight-constituent  tidal  forcing  and 
topographic  wave  drag 

Full  global  3D  output  saved  hourly,  for  one 
month 

Ki.  Oi,  Pi,  and  Qi  (the  four  largest  diurnal  constituents),  as  well 
as  the  same  wind-  and  buoyancy-forcing  used  in  9.7  (Control). 

In  many  applications  involving  tides,  it  is  desirable  to  store 
information  hourly.  However,  at  the  vertical  and  horizontal  resolu¬ 
tions  utilized  here,  it  is  impractical  to  save  full  global  three-dimen¬ 
sional  output  hourly  for  all  five  years  of  the  14.1  simulation.  We 
did  save  daily  25-h  averages  of  three-dimensional  14.1  (All tides) 
output.  For  the  full  five-year  duration  of  14.1  (All tides),  we  saved 
global  hourly  output  of  SSH,  and  other  surface  fields.  We  also  saved 
hourly  full  three-dimensional  output,  over  the  last  three  years  of 
the  run,  in  a  few  domains  of  great  interest  for  the  study  of  internal 
tides,  such  as  Hawai’i,  the  Indonesian  Archipelago,  and  others.  Fi¬ 
nally,  simulation  14.2  (“3Doutput")  is  a  twin  of  14.1  (Alltides)  for 
May  2004  but  saves  full  three-dimensional  model  output  hourly, 
over  the  entire  globe. 

The  combined  size  of  the  stored  output  of  HYCOM  1 4. 1  ( Al  Itides) 
is  68  terabytes  (TB).  and  the  single  month  in  14.2  (3Doutput)  ac¬ 
counts  for  another  3  TB.  This  is  an  enormous  amount  of  material 
to  analyze,  and  we  have  just  begun  to  go  through  our  results.  Thus 
far  a  harmonic  analysis  of  14.1  (Alltides),  commonly  used  to  sepa¬ 
rate  the  contributions  of  the  various  tidal  constituents,  has  been 
performed  at  only  a  limited  number  of  locations,  the  102  pelagic 
tide  gauges  of  Shum  et  al.  (1997).  The  harmonic  analysis  is  used 
to  determine  the  RMS  surface  elevation  errors  of  the  eight  constit¬ 
uents  in  14.1  (Alltides)  with  respect  to  the  tide  gauge  data.  Har¬ 
monic  analysis  on  every  gridpoint  in  such  a  large  model  is  a  very 
time-consuming  endeavor.  For  this  reason,  we  defer  some  of  the 
analyses  we  wish  to  pursue  on  14.1  (Alltides)  to  later  papers.  In  this 
paper,  we  will  show  (1)  results  from  the  harmonic  analysis  of  14.1 
(Alltides)  at  the  tide-gauge  stations,  (2)  other  results  from  14.1  (All¬ 
tides)  which  do  not  require  a  time-consuming  harmonic  analysis, 
and  (3)  some  results  from  harmonic  analysis  of  M2  in  one  day  of 
output  (taken  from  the  second  month,  to  allow  for  spinup)  from 
experiment  14.0  (M2tide).  These  latter  results  are  possible  because 
14.0  (M2tide)  does  not  contain  any  other  tidal  constituents.  How¬ 
ever,  they  should  be  regarded  as  preliminary  because  the  internal 
tide  is  not  necessarily  stationary,  so  that  one  day  of  output  may 
not  be  sufficient  for  a  rigorous  analysis  of  the  internal  tides. 

5.2.  RMS  surface  elevation  errors 

Table  3  shows  the  RMS  time-averaged  elevation  signals  of  the 
eight  largest  consituents  averaged  over  the  102  pelagic  tide  gauges, 
the  elevation  errors  of  year  2004  from  HYCOM  14.1  (Alltides)  with 
respect  to  the  tide  gauge  records  of  these  eight  constituents,  and 
the  percent  of  the  tide  gauge  sea  surface  elevation  variance  of  these 
constituents  captured  by  the  model.  ACHS  may  be  consulted  for 
details  on  how  such  calculations  are  performed.  We  also  analyzed 
years  2003  (last  152  days)  and  2006,  and  came  up  with  virtually 

Table  3 

Time  and  station-averaged  sea-surface  height  (SSH)  signals  at  the  set  of  102  pelagic 
tide  gauges  used  in  Shum  et  al.  (1997X  and  sea-surface  elevation  errors  of  year  2004 
of  our  HYCOM  14.1  eight-constituent  forward  simulation  (Alltides)  with  respect  to 
the  gauges.  Numbers  in  parentheses  denote  percentage  of  SSH  variance  at  the  gauges 
captured  by  HYCOM.  RSS  denotes  root-sum-square  computed  over  all  eight 
constituents. 


Constituent 

Signal  (cm) 

HYCOM  14.1  error  (cm) 

Qi 

1.62 

0.68(82.1) 

0, 

7.76 

2.48  (89.7) 

Pi 

3.62 

0.79  (95.2) 

K, 

11.26 

2.48  (95.1) 

N2 

6.86 

1.40(95.9) 

M2 

33.22 

8.26  (93.8) 

S2 

12.62 

5.17(832) 

3.43 

1.65  (76.9) 

RSS 

39.04 

10.63  (92.6) 
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identical  elevation  errors.  The  overall  percent  variance  captured, 
92.6%,  is  slightly  lower  than  that  captured  in  the  optimally  tuned 
two-layer  simulations  of  ACHS.  The  higher  horizontal  resolution 
used  here  should  improve  the  solutions  (Egbert  et  al.,  2004;  Arbic 
et  al.,  2008).  However,  in  the  latter  models  the  full  spherical  har¬ 
monic  computation  of  SAL  was  utilized,  whereas  here  we  have 
used  only  the  scalar  approximation.  We  conclude  that  for  our  first 
attempt  at  a  mixed  wind-plus-tides  simulation  the  errors  are  rea¬ 
sonably  small.  Based  on  the  experience  in  the  literature  we  believe 
these  errors  will  reduce  with  a  more  rigorous  treatment  of  SAL,  and 
with  the  introduction  of  data  assimilation. 

5.3.  Near-bottom  speeds  and  SSH  voriobility  of  non-tidol  morions 

Because  our  parameterized  topographic  wave  drag  scheme  acts 
on  near-bottom  flows,  and  because  the  25-h  filter  we  utilize  along 
with  the  wave  drag  is  an  imperfect  discriminator  of  tidal  versus 
non-tidal  flows,  it  is  important  to  check  that  non-tidal  near-bottom 
flows  are  not  strongly  affected  with  the  addition  of  tides  and  topo¬ 
graphic  wave  drag.  Fig.  2a  is  a  map  of  the  mean  kinetic  energy  50  m 
above  the  bottom  for  HYCOM  9.7  (Control),  averaged  over  2006. 
Fig.  2b  displays  the  mean  kinetic  energy  50  m  above  the  bottom 
for  non-tidal  flows  in  14.1  (i.e.  based  on  25-h  averaged  daily  cur¬ 


rents  from  Alltides),  also  for  2006.  The  two  figures  were  not  com¬ 
puted  in  exactly  the  same  way,  since  in  Fig.  2a  the  non-tidal  flows 
(i.e.  the  total  flows)  were  saved  as  daily  snapshots  whereas  in 
Fig.  2b  the  non-tidal  flows  were  saved  as  25-h  averages.  However, 
as  shown  in  Arbic  et  al.  (2009),  in  present-day  high-resolution 
models  the  non-tidal  flows  seem  to  be  relatively  unaffected  by 
subsampling  on  scales  of  about  a  day.  Comparison  of  the  two  fig¬ 
ures,  both  of  which  mask  out  regions  shallower  than  1 000  m.  dem¬ 
onstrates  that  on  the  whole,  adding  tides  and  parameterized 
topographic  wave  drag  to  the  model  does  not  reduce  the  non-tidal 
near-bottom  flow.  Indeed,  in  many  regions  it  appears  that  the  non- 
tidal  motions  are  on  the  contrary  stronger  in  the  tidal ly-forced  case 
with  topographic  drag  (14.1;  Alltides)  than  they  are  in  the  non-ti¬ 
dal  case  (9.7;  Control).  We  speculate  that  this  may  be  because,  as 
noted  in  Section  1.  in  the  tidal  simulation  quadratic  bottom  bound¬ 
ary  layer  drag  in  the  deep  ocean  is  effectively  weaker  than  it  is  in 
the  non-tidal  simulation.  Another  possibility  is  that  some  contribu¬ 
tion  from  tidal  flow  is  still  present  in  the  25-h  averages  from  14.1 
(Alltides),  which  would  tend  to  increase  mean  kinetic  energy. 
However,  a  test  of  the  eight-constituent  ACHS  run  shows  that  this 
effect  is  likely  to  contribute  a  maximum  of  7  cm^  s*^  to  the  low- 
frequency  kinetic  energy  in  waters  deeper  than  1000  m.  Larger 
differences  than  this  are  seen  in  the  comparison  of  Fig.  2a  and  b. 


FIfr  2-  Annual  mean  kinetic  energy  (cm^/s^)  SO  m  above  the  bottom  from  (a)  daily  snapshots  from  2006  in  HYCOM  experiment  9.7  (ControIX  which  does  not  have  tidal 
forcing,  (b)  daily  25-h  averages  from  2006  in  HYCOM  experiment  14.1  (Alltides),  in  which  forring  of  the  eight  largest  tidal  constituents  is  included.  The  25-h  averages  filter 
out  most  of  the  tidal  component  of  the  near  bottom  velocities.  Regions  shallower  than  1000  m  are  in  grey.  (For  interpretation  of  the  references  to  colour  in  this  figure  legend, 
the  reader  is  referred  to  the  web  version  of  this  article.) 
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Some  of  the  regions  of  large  mean  kinetic  energy  are  also  where 
the  topographic  wave  drag  is  very  strong  (e.g.  the  Gulf  of  Mexico 
and  some  of  the  Indonesian  Seas).  Our  25-h  averaging  scheme  for 
applying  this  drag  only  to  the  tidal  component  appears  to  induce 
artificially  large  mean  velocities  in  some  locations.  We  are  explor¬ 
ing  alternative  approaches  for  implementing  topographic  wave 
drag  in  subsequent  simulations  with  tides  and  eddies. 

Another  check  on  the  impact  of  including  tides  and  parame¬ 
terized  wave  drag  on  low-frequency  flows  is  given  in  Fig.  3. 
Fig.  3  shows  the  RMS  low-frequency  SSH  variability  over  2004- 
2007  from  HYCOM  9.7  (Control)  and  HYCOM  14.1  (Alltides). 
The  fields  from  the  two  simulations  are  remarkably  similar,  indi¬ 
cating  that  the  low-frequency  SSH  variability  is  not  greatly  af¬ 
fected  by  the  addition  of  tides  and  topographic  wave  drag  to 
the  model. 

5.4.  first-order  impact  of  horiiontally  varying  stratification 

Fig.  4a  displays  the  amplitude  of  the  M2  internal  tide  signature 
in  the  steric  SSH  of  HYCOM  experiment  13.1  (’Tideonly”),  which  is 


run  under  conditions  like  those  in  AGHS  and  SHA;  with  a  horizon¬ 
tally  uniform  two-layer  stratification,  and  no  wind-  and  buoyancy¬ 
forcing.  As  in  AGHS  and  SHA,  the  stratification  was  taken  from  typ¬ 
ical  vertical  profiles  in  subtropical  areas,  which  cover  large  areas  of 
the  world  ocean.  However,  these  stratifications  are  very  different 
from  those  in  polar  regions.  As  Padman  et  al.  (2006)  discuss,  the 
internal  tide  activity  along  the  South  Scotia  Ridge  of  the  Southern 
Ocean  in  AGHS  and  SHA  is  almost  certainly  unrealistically  large. 
Fig.  4a  indicates  that  the  internal  tide  activity  in  HYCOM  13.1 
(Tideonly)  along  the  South  Scotia  Ridge  is  also  large.  Large  internal 
tide  activity  can  also  be  seen  in  AGHS.  SHA,  and  HYCOM  1 3.1  (Tide- 
only)  in  other  polar  regions  e.g.  the  Labrador  Sea  and  the  Southern 
Ocean  south  of  Africa.  Fig.  4b  also  displays  the  M2  steric  SSH  ampli¬ 
tude,  but  computed  from  one  day  of  experiment  14.0(M2tide):  the 
wind,  buoyancy-,  and  M2-forced  ’‘warm-up”  experiment.  In  this 
plot  internal  tide  activity  in  the  polar  regions  is  much  weaker,  thus 
demonstrating  a  first-order  effect  of  horizontally  varying  stratifica¬ 
tion  on  the  internal  tide  field.  In  tropical  and  subtropical  regions 
the  internal  tide  activity  is  generally  stronger  in  the  wind-plus- 
tides  simulation  (Fig.  4b)  than  in  the  tide-only  simulation  (Fig.  4a). 
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Fig.  3.  2004-2007  Root  mean  square  (RMS)  low-frequency  sea  surface  height  (SSH)  variability  (cm)  from  (a)  HYCOM  9.7  (Control;  mean  computed  from  daily  snapshots)  and 
(b)  HYCOM  14.1  (Alltides:  mean  computed  from  daily  25-h  averages). 
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Fig.  4.  Amplitude  (cm)  of  internal  tide  signature  In  steric  SSH  of  HYCOM 
experiment  (a)  13.1  (TIdeonly;  two-layer,  horizontally  uniform  stratification.  M2 
forcing  only),  (b)  14.0  (M2tide;  short  warm-up  run  for  14.1:  32-layer,  horizontally 
non-uniform  stratification,  wind-,  buoyancy-,  and  M2-forcing  included.) 


5.5.  Comparison  of  modeled  interna/  tide  to  sate//ite  a/timeter  data 

We  now  compare  the  modeled  sea  surface  signature  of  internal 
tides  in  the  vicinity  of  Hawai’i  to  the  signatures  seen  in  along-track 
satellite  altimeter  data.  The  altimeter  data  was  obtained  by  per¬ 
sonal  communication  with  Richard  Ray  in  2006.  and  is  an  updated 
version  of  the  data  reported  on  by  Ray  and  Mitchum  (1996, 1997). 
Fig.  5  shows  the  altimeter  tracks  used  in  the  comparison.  The  blue 
lines  in  Fig.  6  show  the  M2  elevation  amplitudes  and  phases  along 
track  number  125,  in  observations  and  in  HYCOM  14.0  (M2tide). 
The  red  lines  denote  the  low-pass  filtered  (barotropic)  versions  of 
the  full  (blue  lines)  signal.  In  Fig.  7  we  display  the  difference  be¬ 
tween  the  blue  and  red  lines  in  Fig.  6,  i.e.  the  perturbations  to 
the  M2  elevation  amplitudes  and  phases  at  the  sea  surface  due  to 
internal  tides.  The  peaks  and  troughs  in  the  modeled  perturbations 
clearly  have  similar  amplitude  and  horizontal  length  scale  to  those 
in  the  observed  perturbations,  but  equally  clearly  do  not  match  the 
peaks  and  troughs  in  the  observed  perturbations  along  the  entire 
track  length,  especially  far  from  the  source  at  Hawai’i.  In  contrast, 
when  high-resolution  regional  models  forced  at  their  horizontal 
boundaries  by  altimeter-derived  tidal  amplitudes  are  compared 
to  altimeter  data  (e.g.  Carter  et  al.,  2008  and  references  therein, 
among  several),  the  comparison  is  better.  The  RMS  of  the  internal 
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Fig.  5.  Tracks  for  which  satellite  altimeter  data  (around  Hawaii)  on  surface 
signature  of  M2  internal  tides  is  utilized  here. 


tide  perturbations,  averaged  over  all  of  the  tracks  shown  in  Fig.  5 
(using  the  latitude  and  longitude  bounds  shown  in  the  Figure), 
are  given  in  Table  4.  RMS  values  are  given  for  observations,  AGHS, 
SHA,  and  HYCOM  14.0  {M2tide).  The  RMS  of  the  differences  be¬ 
tween  the  observed  and  modeled  perturbations  (i.e.  differences  be¬ 
tween  the  blue  and  red  curves  in  Fig.  7)  are  given  in  parentheses  in 
Table  4.  The  AGHS  internal  tides  are  too  weak,  probably  because  of 
the  relatively  low  1/4®  resolution  used  there,  as  evidenced  by  the 
low  RMS  values  compared  to  observations.  The  magnitudes  of 
the  SHA  and  HYCOM  14.0  (M2tide)  pertubations  are  closer  to  those 
seen  in  observations,  for  both  amplitude  and  phase.  The  poor 
match  of  the  peaks  and  troughs  in  the  modeled  and  observed  per¬ 
turbations  is  seen  in  the  values  of  the  RMS  differences  between 
perturbations  (the  parenthetical  values),  which  are  nearly  as  large 
or  larger  than  the  RMS  values  seen  in  the  observations. 


(a)  Observed  amplitude 


(c)  Observed  phase 


Fig.  6.  Amplitudes  (a.  observed  and  b.  HYCOM  I4.0-M2tide)  and  phases  (c. 
observed  and  d.  HYCOM  14.0-M2tide)  of  the  Mj  sea  surface  elevation  along 
altimetric  track  number  125.  HYCOM  14.0  (M2tide)  is  a  32-layer,  wind-,  buoyancy-, 
and  M2-  forced  simulation.  Blue  lines  represent  full  signal  (barotropic  plus 
baroclinicX  red  lines  represent  low-pass  filtered  (barotropic)  signal. 
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(a)  Track  125  amplitude  perturbation 


(b)  Track  125  phase  perturbation 


Fig.  7.  Internal  tide  perturbations  to  (a)  amplitude  and  (b)  phase  of  Mj  sea  surface 
elevation  along  altimetric  track  125.  Altimetric  observations  are  in  blue,  while 
HYCOM  14.0  (M2tide)  is  in  red.  (For  interpretation  of  the  references  to  colour  in  this 
figure  legend,  the  reader  is  referred  to  the  web  version  of  this  article.) 


Table  4 

RMS  of  the  internal  tide  perturbations  to  M2  sea  surface  elevation  amplitudes  and 
phases,  computed  across  all  of  the  tracks  shown  in  Fig.  5.  from  altimetric 
observations,  AGHS,  SHA,  and  HYCOM  14.0  (M2tide).  RMS  values  of  the  difference 
in  perturbations  (model  minus  observations)  are  given  in  parentheses. 


Source 

RMS  amplitude 
perturbation  (cm) 

RMS  phase  perturbation 
(degrees) 

Observations 

0.87 

4.35 

AGHS 

0.40  (0.86) 

1.91  (3.93) 

SHA 

1.07(1.29) 

4.66  (5.64) 

HYCOM  14.0(M2tide) 

1.03  (1.15) 

4.42  (4.58) 

5.6.  Co-existence  of  tides  and  the  eddying  general  circulation 

We  now  show  some  figures  which  visually  demonstrate  the  co¬ 
existence  of  tides  and  the  eddying  general  circulation  in  the  HY¬ 
COM  simulations.  In  Fig.  8  we  show  snapshots  of  the  non-steric 
and  steric  sea-surface  heights  in  the  Southwest  Padfic  sector  of 
HYCOM  14.1  (Alltides).  The  non-steric  height  (Fig.  8a)  is  dominated 
by  the  large-scale  barotropic  tide.  The  eddying  general  circulation 
features  smaller  horizontal  scales  and  is  easily  discernible  in  Fig.  8b 
(the  steric  plot),  as  in  many  previous  studies  of  high-resolution 
ocean  models.  Internal  tides  are  visible  as  small-scale  speckled  pat¬ 
terns  in  the  upper  portion  of  the  steric  plot.  It  is  difficult  by  eye  to 


Fig.  8.  Snapshot  of  (a)  non-steric  and  (b)  steric  sea  surface  heights  (m)  in  the  Southwest  Pacific  on  June  30.  2006  at  OOZ,  from  HYCOM  14.1  (Alltidcs). 
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Fig.  9.  Global  difference  in  steric  sea  surface  heights  (m)  from  snapshots  taken  6  h 
apart;  June  30.  2006.  06Z-00Z,  from  HYCOM  14.1  (Alltides). 

discern  differences  in  steric  height  fields  lying  close  together  in 
time.  A  global  map  of  differences  between  two  steric  height  fields 
6  h  apart  is  displayed  in  Fig.  9.  Even  with  a  color  scale  that  covers  a 
range  5x  smaller,  the  meso-  and  gyre-scale  general  circulation  fea¬ 
tures  in  Fig.  8b  are  absent  in  the  difference  plot.  Instead,  we  see  the 


much  higher  frequency  internal  tides,  which  show  up  as  beams  as 
in  Fig.  4. 

The  co-existence  of  tides  and  the  eddying  general  circulation 
can  be  seen  more  easily  in  animations  which  we  have  submitted 
along  with  this  paper,  Hawaii.fli  (Supplementary  mmcl.avi)  shows 
the  steric  and  non-steric  sea  surface  heights  in  a  region  around  Ha¬ 
waii,  for  the  last  five  days  of  June  2004.  The  non-steric  field,  dom¬ 
inated  by  the  barotropic  tide,  evolves  rapidly  in  time  and  features 
large  horizontal  scales.  In  the  steric  field,  which  features  smaller 
horizontal  scales,  the  higher-frequency  internal  tide  signals  course 
rapidly  through  the  eddying  general  circulation,  which  appears  to 
be  at  a  standstill  on  these  short  timescales.  Stericssh.gif  (Supple¬ 
mentary  mmc2.gif).  a  movie  of  the  steric  SSH  covering  a  much  lar¬ 
ger  area  as  well  as  a  longer  time  period,  shows  that  internal  tides 
are  ubiquitous  throughout  the  world  ocean. 

Finally,  we  give  some  indication  of  the  vertical  structure  of  the 
simulation  in  Fig.  10a  and  b.  These  figures,  computed  from  exper¬ 
iment  14.1  (Alltides),  display  the  zonal  component  of  velocity  (u) 
in  the  upper  waters  of  a  meridional  section  running  through  Ha¬ 
waii.  Fig,  10b  shows  the  25-h  mean  while  Fig.  10a  shows  the  snap¬ 
shot  at  noon  Zulu  (UST).  Layer  interfaces  are  shown  as  solid  black 
lines,  and  the  thick  black  line  is  the  mixed  layer  depth.  The  hybrid 
nature  of  HYCOM’s  vertical  coordinate  is  illustrated  by  the  increas¬ 
ing  number  of  near-surface  layers  that  are  flat  (i.e.  in  z-coordi- 
nates)  the  further  north  in  the  plot.  However,  the  majority  of  the 
layers  are  isopycnal  and  so  give  an  indication  of  the  density  struc¬ 
ture.  There  is  much  more  structure  in  the  small-scale  perturbations 
of  the  interfaces  between  isopycnal  layers  in  Fig.  10a  than  in 
Fig.  1  Ob,  indicating  that  many  of  the  perturbations  are  due  to  inter¬ 
nal  tides.  These  perturbations  involve  changes  in  vertical  density 


Fig.  10.  Zonal  component  of  velocity  (u;  cm  s  ')  in  ISS^W  section  through  Hawai’i  on  June  30. 2006  from  HYCOM  14.1  (Alltides):  (a)  snapshot  and  (b)  25-h  mean.  Isopycnal 
locations  (black  lines)  shown  versus  depth  in  meters  (right  axes). 
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Structure  and  will  have  signatures  in  the  steric  SSH.  Likewise,  there 
is  much  more  vertical  structure  in  the  velocity  field  in  the  snapshot 
than  in  the  25-h  average,  indicating  that  the  tides  are  a  strong  sig¬ 
nal.  and  have  significant  vertical  structure,  in  that  field  as  well. 

6.  Summary  and  discussion 

In  this  paper  we  have  shown  some  preliminary  results  of  HY- 
COM  simulations  which  simultaneously  resolve  barotropic  tides, 
baroclinic  tides,  and  an  eddying  general  circulation.  The  nominal 
horizontal  resolution  of  the  simulation  is  1/12.5°,  and  there  are 
32  hybrid  layers  in  the  vertical  direction.  We  have  shown  that  a 
parameterized  topographic  wave  drag  can  be  inserted  which  yields 
a  reasonably  accurate  sea  surface  elevation  of  the  barotropic  tide  at 
the  same  time  that  the  near-bottom  flows  of  non-tidal  motions, 
and  low-frequency  SSH  variability,  are  not  strongly  affected.  The 
barotropic  tide  in  the  baroclinic  simulations  presented  here  is  of 
comparable  accuracy  to  that  in  the  main  baroclinic  simulations 
of  Arbic  et  al.  (2004  -  AGHS),  and  is  considerably  more  accurate 
than  that  in  the  main  baroclinic  simulation  of  Simmons  et  al. 
(2004  -  SHA). 

The  stratification  in  the  simulations  presented  here  can  vary  in 
the  horizontal  direction,  since  wind-  and  buoyancy  forcing  is  pres¬ 
ent  to  support  such  variations.  In  contrast,  the  stratification  in  the 
earlier  global  baroclinic  tide  simulations  of  AGHS  and  SHA,  which 
did  not  include  wind-  and  buoyancy-forcing,  was  horizontally  uni¬ 
form.  In  AGHS  and  SHA  a  typical  midlatitude  stratification  was 
used  throughout  the  entire  globe,  and  internal  tide  activity  in  some 
polar  regions  (for  instance,  the  Scotia  Sea)  was  almost  certainly 
artificially  high  (Padman  et  al..  2006).  Comparison  of  the  internal 
tide  signature  at  the  sea  surface  in  HYCOM  runs  with  a  horizontally 
uniform  stratification  and  tidal  forcing  only  versus  the  more  realis¬ 
tic  horizontally  varying  stratification  in  a  wind,  buoyancy,  and  tid- 
ally  forced  run.  indicates  that  internal  tide  activity  in  polar  regions 
is  much  reduced  in  the  latter  compared  to  the  former.  Thus  the 
allowance  of  a  horizontally  varying  stratification  with  the  inclusion 
of  wind-  and  buoyancy-forcing  has  a  first-order  effect  on  the  inter¬ 
nal  tide  field. 

Preliminary  comparisons  of  the  surface  signature  of  the  Ma 
internal  tide  in  the  region  around  Hawai'i  with  satellite  altimeter 
data  indicate  that  the  internal  tides  in  HYCOM  appear  to  have 
approximately  correct  magnitude.  Similar  comparisons  show  that 
the  AGHS  internal  tides  are  too  weak,  while  the  SHA  internal  tides 
are  of  similar  amplitude  to  the  HYCOM  internal  tides.  Thus  the  HY¬ 
COM  simulations  presented  here,  unlike  either  the  AGHS  or  SHA 
simulations,  contain  reasonably  accurate  barotropic  tides,  and 
baroclinic  tides  that  are  at  least  of  the  correct  magnitude,  at  the 
same  time.  However,  there  is  not  in  general  a  close  match  between 
the  peaks  and  troughs  seen  in  the  internal  tide  perturbations  to  sea 
surface  elevation  in  the  observations  and  in  the  HYCOM  model.  In  a 
planned  future  paper  we  will  investigate  the  surface  signature  of 
the  internal  tides,  and  their  comparison  to  satellite  altimeter  data, 
in  much  more  detail. 
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